Comparison of rRNA depletion methods for efficient bacterial mRNA sequencing.

Current methods of high-throughput RNA sequencing of prokaryotes, including transcriptome analysis or ribosomal profiling, need deep sequencing to achieve sufficient numbers of effective reads (e.g., mapping to mRNA) in order to also find weakly expressed genetic elements. The fraction of high-quality reads mapping to coding RNAs (i.e., mRNA) is mainly influenced by the large content of rRNA and, to a lesser extent, tRNA in total RNA. Thus, depletion of rRNA increases coverage and thus sequencing costs. RiboZero, a depletion kit based on probe hybridisation and rRNA-removal was found to be most efficient in the past, but it was discontinued in 2018. To facilitate comparability with previous experiments and to help choose adequate replacements, we compare three commercially available rRNA depletion kits also based on hybridization and magnetic beads, i.e., riboPOOLs, RiboMinus and MICROBExpress, with the former RiboZero. Additionally, we constructed biotinylated probes for magnetic bead capture and rRNA depletion in this study. Based on E. coli, we found similar efficiencies in rRNA depletion for riboPOOLs and the self-made depletion method; both comparable to the former RiboZero, followed by RiboMinus, succeeded by MICROBExpress. Further, our in-house protocol allows customized species-specific rRNA or even tRNA depletion or depletion of other RNA targets. Both, the self-made biotinylated probes and riboPOOLs, were most successful in reducing the rRNA content and thereby increasing sequencing depth concerning mRNA reads. Additionally, the number of reads matching to weakly expressed genes are increased. In conclusion, the self-made specific biotinylated probes and riboPOOLs are an adequate replacement for the former RiboZero. Both are very efficient in depleting rRNAs, increasing mRNA reads and thus sequencing efficiency.

Facklamia lactis sp. nov., isolated from raw milk.

Two strains of a Gram-staining-positive species were isolated from German bulk tank milk. On the basis of their 16S rRNA sequences they were affiliated to the genus Facklamia but could not be assigned to any species with a validly published name. Facklamia miroungae ATCC BAA-466T (97.3 % 16S rRNA sequence similarity), Facklamia languida CCUG 37842T (96.9 %), and Facklamia hominis CCUG 36813T (96.6 %) are the closest relatives. In the 16S rRNA phylogeny and in the core-genome phylogeny strains WS 5301T and WS 5302 form a well-supported, separate lineage. Pairwise average nucleotide identity calculated using MUMmer (ANIm) between WS 5301T and type strains of other Facklamia species is well below the species cut-off (95 %) and ranges from 83.4 to 87.7 %. The DNA G+C content of the type strain is 36.4 mol% and the assembly size of the genome is 2.2 Mb. Cells of WS 5301T are non-motile, non-endospore-forming, oxidase-negative, catalase-negative and facultatively anaerobic cocci. The fastidious species grows at 10-40 °C and with up to 7.0 % (w/v) NaCl in BHI supplemented with 5 g l-1 yeast extract. Major polar lipids are phosphatidylglycerol, diphosphatidylglycerol and two glycolipids. Predominant fatty acids are C16 : 1ω9c and C18 : 1ω9c. On the basis of their genomic, physiological and chemotaxonomic characteristics the strains examined in this study represent the same, hitherto unknown species. We propose the name Facklamia lactis sp. nov. for which WS 5301T (=DSM 111018T=LMG 31861T) is the type strain and WS 5302 (=DSM 111019=LMG 31862) is an additional strain of this novel species.

Simultaneous quantification of the most common and proteolytic Pseudomonas species in raw milk by multiplex qPCR.

The heat-stable peptidase AprX, secreted by psychrotolerant Pseudomonas species in raw milk, is a major cause of destabilization and premature spoilage of ultra-high temperature (UHT) milk and milk products. To enable rapid detection and quantification of seven frequent and proteolytic Pseudomonas species (P. proteolytica, P. gessardii, P. lactis, P. fluorescens, P. protegens, P. lundensis, and P. fragi) in raw milk, we developed two triplex qPCR assays taking into account species-dependent differences in AprX activity. Besides five species-specific hydrolysis probes, targeting the aprX gene, a universal rpoB probe was included in the assay to determine the total Pseudomonas counts. For all six probes, linear regression lines between Cq value and target DNA concentration were obtained in singleplex as well as in multiplex approaches, yielding R2 values of > 0.975 and amplification efficiencies of 85-97%. Moreover, high specificity was determined using genomic DNA of 75 Pseudomonas strains, assigned to 57 species, and 40 other bacterial species as templates in the qPCR. Quantification of the target species and total Pseudomonas counts resulted in linear detection ranges of approx. 103-107 cfu/ml, which correspond well to common Pseudomonas counts in raw milk. Application of the assay using 60 raw milk samples from different dairies showed good agreement of total Pseudomonas counts calculated by qPCR with cell counts derived from cultivation. Furthermore, a remarkably high variability regarding the species composition was observed for each milk sample, whereby P. lundensis and P. proteolytica/P. gessardii were the predominant species detected. KEY POINTS: • Multiplex qPCR for quantification of seven proteolytic Pseudomonas species and total Pseudomonas counts in raw milk • High specificity and sensitivity via hydrolysis probes against aprX and rpoB • Rapid method to determine Pseudomonas contamination in raw milk and predict spoilage potential.

Pseudomonas cremoris sp. nov., a novel proteolytic species isolated from cream.

During a study investigating the microbiota of raw milk and its semi-finished products, strains WS 5106T and WS 5096 were isolated from cream and skimmed milk concentrate. They could be assigned to the genus Pseudomonas by their 16S rRNA sequences, but not to any validly named species. In this work, a polyphasic approach was used to characterize the novel strains and to investigate their taxonomic status. Examinations based on the topology of core genome phylogenomy as well as average nucleotide identity (ANIm) comparisons suggested a novel Pseudomonas species within the Pseudomonas fluorescens subgroup. With pairwise ANIm values of 90.1 and 89.8 %, WS 5106T was most closely related to Pseudomonas nabeulensis CECT 9765T and Pseudomonas kairouanensis CECT 9766T. The G+C content of strain WS 5106T was 60.1 mol%. Morphologic analyses revealed Gram-stain-negative, aerobic, catalase and oxidase positive, rod-shaped and motile cells. Proteolysis on skimmed milk agar as well as lipolysis on tributyrin agar occurred at both 28 and 6 °C. Tolerated growth conditions were temperatures between 4 and 34 °C, pH values between 6.0 and 8.0, and salt concentrations of up to 5 %. Fatty acid profiles showed a pattern typical for Pseudomonas, with C16 : 0 as the dominant component. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol and the dominating quinone was Q-9. Based on these results, it is proposed to classify the strains as a novel species, Pseudomonas cremoris sp. nov., with WS 5106T (=DSM 111143T=LMG 31863T) as type strain and WS 5096 (=DSM 111129=LMG 31864) as an additional strain.

High counts of thermophilic spore formers in dairy powders originate from persisting strains in processing lines.

During the last decades, thermophilic spore counts became a very important quality parameter for manufacturers with regard to powdered dairy products. Low-spore count powders are highly demanded but challenging to produce when high production volume and long process times are intended. In this study a detailed monitoring of microbial levels in three skim milk powder plants was conducted. Anoxybacillus flavithermus was found to be primarily responsible for increased spore levels with increasing spore numbers being detected after 6-8 h already during initial processing steps. Simultaneously, the species composition shifted from a diverse bulk tank milk microbiota where different Bacillus species represented around 90% of the thermophilic bacteria to a dominance of A. flavithermus in the end product. The analysis of A. flavithermus isolates from different powder batches with RAPD PCR revealed recurring patterns in each of the eight German manufacturers sampled over several months. The high relatedness of isolates exhibiting identical RAPD patterns was exemplified by cgMLST based on whole genome sequences. We assume that A. flavithermus strains persisted in production plants and were not eliminated by cleaning. It is concluded that such persisting strains recurrently recontaminated subsequent powder productions. The data highlight that a targeted optimization of cleaning and disinfection procedures is the most promising measure to effectively reduce thermophilic spore counts in German dairy powders.

Fundicoccus ignavus gen. nov., sp. nov., a novel genus of the family Aerococcaceae isolated from bulk tank milk.

Three strains of a Gram-stain-positive, catalase-negative, facultative anaerobic, and coccoid species were isolated from German bulk tank milk. Phylogenetic analyses based on the 16S rRNA gene sequences indicated that the three strains (WS4937T, WS4759 and WS5303) constitute an independent phylogenetic lineage within the family Aerococcaceae with Facklamia hominis CCUG 36813T (93.7-94.1 %) and Eremococcus coleocola M1831/95/2T (93.5 %) as most closely related type species. The unclassified strains demonstrated variable growth with 6.5 % (w/v) NaCl and tolerated pH 6.5-9.5. Growth was observed from 12 to 39 °C. Their cell-wall peptidoglycan belongs to the A1α type (l-Lys-direct) consisting of alanine, glutamic acid and lysine. The predominant fatty acids were C16 : 1 ω9c, C16 : 0 and C18 : 1 ω9c and in the polar lipids profile three glycolipids, a phospholipid, phosphatidylglycerol, phosphoglycolipid and diphosphatidylglycerol were found. The G+C content of strain WS4937T was 37.4 mol% with a genome size of ~3.0 Mb. Based on phylogenetic, phylogenomic and biochemical characterizations, the isolates can be demarcated from all other genera of the family Aerococcaceae and, therefore, the novel genus Fundicoccus gen. nov. is proposed. The type species of the novel genus is Fundicoccus ignavus gen. nov., sp. nov. WS4937T (=DSM 109652T=LMG 31441T).

Genetic Organization of the aprX-lipA2 Operon Affects the Proteolytic Potential of Pseudomonas Species in Milk.

Psychrotolerant Pseudomonas species are a main cause of proteolytic spoilage of ultra-high temperature (UHT) milk products due to the secretion of the heat-resistant metallopeptidase AprX, which is encoded by the first gene of the aprX-lipA2 operon. While the proteolytic property has been characterized for many different Pseudomonas isolates, the underlying aprX-lipA2 gene organization was only described for a few strains so far. In this study, the phylogenomic analysis of 185 Pseudomonas type strains revealed that the presence of aprX is strongly associated to a monophylum composed of 81 species, of which 83% carried the aprX locus. Furthermore, almost all type strains of known milk-relevant species were shown to be members of the three monophyletic groups P. fluorescens, P. gessardii, and P. fragi. In total, 22 different types of aprX-lipA2 genetic organizations were identified in the genus, whereby 31% of the species tested carried the type 1 operon structure consisting of eight genes (aprXIDEF prtAB lipA2). Other genetic structures differed from type 1 mainly in the presence and location of genes coding for two lipases (lipA1 and lipA2) and putative autotransporters (prtA and prtB). The peptidase activity of 129 strains, as determined on skim milk agar and in UHT-milk, correlated largely with different aprX-lipA2 gene compositions. Particularly, isolates harboring the type 1 operon were highly proteolytic, while strains with other operon types, especially ones lacking prtA and prtB, exhibited significantly lower peptidase activities. In conclusion, the phylogenomic position and the aprX-lipA2 gene organization specify the proteolytic potential of Pseudomonas isolates. In addition, however, an interplay of several environmental factors and intrinsic traits influences production and activity of AprX, leading to strain-specific proteolytic phenotypes.

Recommendations for bacterial ribosome profiling experiments based on bioinformatic evaluation of published data.

Ribosome profiling (RIBO-Seq) has improved our understanding of bacterial translation, including finding many unannotated genes. However, protocols for RIBO-Seq and corresponding data analysis are not yet standardized. Here, we analyzed 48 RIBO-Seq samples from nine studies of Escherichia coli K12 grown in lysogeny broth medium and particularly focused on the size-selection step. We show that for conventional expression analysis, a size range between 22 and 30 nucleotides is sufficient to obtain protein-coding fragments, which has the advantage of removing many unwanted rRNA and tRNA reads. More specific analyses may require longer reads and a corresponding improvement in rRNA/tRNA depletion. There is no consensus about the appropriate sequencing depth for RIBO-Seq experiments in prokaryotes, and studies vary significantly in total read number. Our analysis suggests that 20 million reads that are not mapping to rRNA/tRNA are required for global detection of translated annotated genes. We also highlight the influence of drug-induced ribosome stalling, which causes bias at translation start sites. The resulting accumulation of reads at the start site may be especially useful for detecting weakly expressed genes. As different methods suit different questions, it may not be possible to produce a "one-size-fits-all" ribosome profiling data set. Therefore, experiments should be carefully designed in light of the scientific questions of interest. We propose some basic characteristics that should be reported with any new RIBO-Seq data sets. Careful attention to the factors discussed should improve prokaryotic gene detection and the comparability of ribosome profiling data sets.

Brevilactibacter flavus gen. nov., sp. nov., a novel bacterium of the family Propionibacteriaceae isolated from raw milk and dairy products and reclassification of Propioniciclava sinopodophylli as Brevilactibacter sinopodophylli comb. nov.

Eight facultatively anaerobic rod-shaped bacteria were isolated from raw milk and two other dairy products. Results of phylogenetic analyses based on 16S rRNA gene sequences showed that the isolates are placed in a distinct lineage within the family Propionibacteriaceae with Propioniciclava sinopodophylli and Propioniciclava tarda as the closest relatives (94.6 and 93.5 % similarity, respectively). The cell-wall peptidoglycan contained meso-diaminopimelic acid, alanine and glutamic acid and was of the A1γ type (meso-DAP-direct). The major cellular fatty acid was anteiso-C15 : 0 and the major polar lipids were diphosphatidylglycerol, phosphatidyglycerol and three unidentified glycolipids. The quinone system contained predominantly menaquinone MK-9(H4). The G+C content of the genomic DNA of strain VG341T was 67.7 mol%. The whole-cell sugar pattern contained ribose, rhamnose, arabinose and galactose. On the basis of phenotypic and genetic data, eight strains (VG341T, WS4684, WS4769, WS 4882, WS4883, WS4901, WS4902 and WS4904) are proposed to be classified as members of a novel species in a new genus of the family Propionibacteriaceae, for which the name Brevilactibacter flavus gen. nov., sp. nov. is proposed. The type strain is VG341T (=WS4900T=DSM 100885T=LMG 29089T) and seven additional strains are WS4684, WS4769, WS4882, WS4883, WS4901, WS4902 and WS4904. Furthermore, we propose the reclassification of P. sinopodophylli as Brevilactibacter sinopodophylli comb. nov.

Pseudomonas haemolytica sp. nov., isolated from raw milk and skimmed milk concentrate.

Two strains, WS 5063T and WS 5067, isolated from raw cow's milk and skimmed milk concentrate, could be affiliated as members of the same, hitherto unknown, Pseudomonas species by 16S rRNA and rpoD gene sequences. Multilocus sequence and average nucleotide identity (ANIm) analyses based on draft genome sequences confirmed the discovery of a novel Pseudomonas species. It was most closely related to Pseudomonas synxantha DSM 18928T with an ANIm of 91.4 %. The DNA G+C content of WS 5063T was 60.0 mol %. Phenotypic characterizations showed that the isolates are rod-shaped, motile, catalase- and oxidase-positive, and aerobic. Growth occurred at 4-34 °C and at pH values of pH 5.5-8.0. Both strains showed strong ß-haemolysis on blood agar. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol and diphosphatidylglycerol. The dominant quinone was Q-9 (90 %), but noticeable amounts of Q-8 (9 %) and traces of Q-7 were also detected. Fatty acid profiles were typical for Pseudomonas species and exhibited C16 : 0 as a major component. Based on these results, we conclude that both strains belong to a novel species, for which the name Pseudomonas haemolytica sp. nov. is proposed. The type strain is WS 5063T (=DSM 108987T=LMG 31232T) and an additional strain is WS 5067 (=DSM 108988=LMG 31233).

Pseudomonas saxonica sp. nov., isolated from raw milk and skimmed milk concentrate.

A polyphasic approach was used to investigate the taxonomic status of two bacterial strains, WS 5072T and WS 5092, isolated from skimmed milk concentrate and raw cow's milk. The 16S rRNA and rpoD gene sequences affiliated the strains to the same, hitherto unknown, Pseudomonas species. Further examinations of the draft genomes based on multilocus sequence analysis and average nucleotide identity confirmed the presence of a novel Pseudomonas species. It was most closely related to Pseudomonas fragi DSM 3456T with 86.3 % ANIm. The DNA G+C content of strain WS 5072T was 56.3 mol%. Cells were aerobic, Gram-negative, catalase and oxidase positive, rod-shaped and motile. Growth occurred at 4-34 °C, pH 5.5-8.0 and with salt concentrations of up to 7 %. The major cellular polar lipids were phosphatidylethanolamine, phosphatidylglycerol, and diphosphatidylglycerol. The dominating quinone was Q-9 with 94 %, with noticeable amounts of Q-8 (5 %) and traces of Q-7 and Q-10. Fatty acid profiles showed a composition common for Pseudomonas with the major component C16 : 0. Based on these results, the novel species Pseudomonas saxonica sp. nov. is proposed, with the type strain WS 5072T (=DSM 108989T=LMG 31234T) and the additional strain WS 5092 (=DSM 108990=LMG 31235).

Proposal of Lysobacter pythonis sp. nov. isolated from royal pythons (Python regius).

The bacterial strains 4284/11T and 812/17 isolated from the respiratory tract of two royal pythons in 2011 and 2017, respectively were subjected to taxonomic characterization. The 16S rRNA gene sequences of the two strains were identical and showed highest sequence similarities to Lysobacter tolerans UM1T (97.2%) and Luteimonas aestuarii DSM 19680T (96.7 %). The two strains were identical in the sequences of the 16S-23S rRNA internal transcribed spacer (ITS) and partial groEL gene sequences and almost identical in genomic fingerprints. In the ITS sequence Ly. tolerans DSM 28473T and in the groEL nucleotide sequence Luteimonas mephitis DSM 12574T showed the highest similarity. In silico DDH analyses using genome sequence based ANIb and gANI similarity coefficients demonstrated that strain 4284/11T represents a novel species and revealed Ly. tolerans UM1T as the next relative (ANIb = 76.2 %, gANI = 78.0 %). Based on the topology of a core gene phylogeny strain 4284/11T could be assigned to the genus Lysobacter. Chemotaxonomic characteristics including polyamine pattern, quinone system, polar lipid profile and fatty acid profile were in accordance with the characteristics of the genera Lysobacter and Luteimonas. Strains 4284/11T and 812/17 could be differentiated from the type strains of the most closely related species by several physiological tests. In conclusion we are here proposing the novel species Lysobacter pythonis sp. nov. The type strain is 4284/11T (= CCM 8829T = CCUG 72164T = LMG 30630T) and strain 812/17 (CCM 8830) is a second strain of this species.

Complete Circular Genome Sequence and Temperature Independent Adaptation to Anaerobiosis of Listeria weihenstephanensis DSM 24698.

The aim of this study was to analyze the adaptation of the environmental Listeria weihenstephanensis DSM 24698 to anaerobiosis. The complete circular genome sequence of this species is reported and the adaptation of L. weihenstephanensis DSM 24698 to oxygen availability was investigated by global transcriptional analyses via RNAseq at 18 and 34°C. A list of operons was created based on the transcriptional data. Forty-two genes were upregulated anaerobically and 62 genes were downregulated anaerobically. The oxygen dependent gene expression of selected genes was further validated via qPCR. Many of the differentially regulated genes encode metabolic enzymes indicating broad metabolic adaptations with respect to oxygen availability. Genes showing the strongest oxygen-dependent adaption encoded nitrate (narGHJI) and nitrite (nirBD) reductases. Together with the observation that nitrate supported anaerobic growth, these data indicate that L. weihenstephanensis DSM 24698 performs anaerobic nitrate respiration. The wide overlap between the oxygen-dependent transcriptional regulation at 18 and 34°C suggest that temperature does not play a key role in the oxygen-dependent transcriptional regulation of L. weihenstephanensis DSM 24698.

Pseudomonas lactis sp. nov. and Pseudomonas paralactis sp. nov., isolated from bovine raw milk.

Five strains, designated WS 4672T, WS 4998, WS 4992T, WS 4997 and WS 5000, isolated from bovine raw milk formed two individual groups in a phylogenetic analysis. The most similar species on the basis of 16S rRNA gene sequences were Pseudomonas azotoformans IAM 1603T, Pseudomonas gessardii CIP 105469T and Pseudomonas libanensis CIP 105460T showing 99.7-99.6 % similarity. Using rpoD gene sequences Pseudomonas veronii LMG 17761T (93.3 %) was most closely related to strain WS 4672T and Pseudomonas libanensis CIP 105460T to strain WS 4992T (93.3 %). The five strains could be differentiated from their closest relatives and from each other by phenotypic and chemotaxonomic characterization and ANIb values calculated from draft genome assemblies. ANIb values of strains WS 4992T and WS4671T to the closest relatives are lower than 90 %. The major cellular polar lipids of both strains are phosphatidylethanolamine, phosphatidylglycerol, a phospholipid and diphosphatidylglycerol, and their major quinone is Q-9. The DNA G+C content of strains WS 4992T and WS 4672T were 60.0  and 59.7  mol%, respectively. Based on these genotypic and phenotypic traits two novel species of the genus Pseudomonas are proposed: Pseudomonas lactis sp. nov. [with type strain WS 4992T (=DSM 29167T=LMG 28435T) and the additional strains WS 4997 and WS 5000], and Pseudomonasparalactis sp. nov. [with type strain WS 4672T (=DSM 29164T=LMG 28439T) and additional strain WS 4998].

Draft Genome Sequence of Lysinibacillus xylanilyticus SR-86.

Lysinibacillus xylanilyticus belongs to the family Bacillaceae and was first described in 2010 with the type strain L. xylanilyticus XDB9. It is able to both degrade xylan and use it as the sole carbon source. Here, we describe the 4.8-Mb genome of the strain L. xylanilyticus SR-86, which was isolated from organic waste.

Draft Genome Sequence of the Xanthan Producer Xanthomonas campestris LMG 8031.

Here, we report the draft genome sequence of Xanthomonas campestris LMG 8031, for which nearly no genetic information is available, despite its good xanthan-producing properties. We performed an Illumina-based sequencing approach of LMG 8031. The genome revealed a 5.0-Mb chromosome having 4,434 coding sequences and a G+C content of 65%.

Optimized Illumina PCR-free library preparation for bacterial whole genome sequencing and analysis of factors influencing de novo assembly.

BACKGROUND: Next-generation sequencing (NGS) technology has paved the way for rapid and cost-efficient de novo sequencing of bacterial genomes. In particular, the introduction of PCR-free library preparation procedures (LPPs) lead to major improvements as PCR bias is largely reduced. However, in order to facilitate the assembly of Illumina paired-end sequence data and to enhance assembly performance, an increase of insert sizes to facilitate the repeat bridging and resolution capabilities of current state of the art assembly tools is needed. In addition, information concerning the relationships between genomic GC content, library insert size and sequencing quality as well as the influence of library insert size, read length and sequencing depth on assembly performance would be helpful to specifically target sequencing projects. RESULTS: Optimized DNA fragmentation settings and fine-tuned resuspension buffer to bead buffer ratios during fragment size selection were integrated in the Illumina TruSeq(®) DNA PCR-free LPP in order to produce sequencing libraries varying in average insert size for bacterial genomes within a range of 35.4-73.0 % GC content. The modified protocol consumes only half of the reagents per sample, thus doubling the number of preparations possible with a kit. Examination of different libraries revealed that sequencing quality decreases with increased genomic GC content and with larger insert sizes. The estimation of assembly performance using assembly metrics like corrected NG50 and NGA50 showed that libraries with larger insert sizes can result in substantial assembly improvements as long as appropriate assembly tools are chosen. However, such improvements seem to be limited to genomes with a low to medium GC content. A positive trend between read length and assembly performance was observed while sequencing depth is less important, provided a minimum coverage is reached. CONCLUSIONS: Based on the optimized protocol developed, sequencing libraries with flexible insert sizes and lower reagent costs can be generated. Furthermore, increased knowledge about the interplay of sequencing quality, insert size, genomic GC content, read length, sequencing depth and the assembler used will help molecular biologists to set up an optimal experimental and analytical framework with respect to Illumina next-generation sequencing of bacterial genomes.

Draft Genome Sequences of Three European Laboratory Derivatives from Enterohemorrhagic Escherichia coli O157:H7 Strain EDL933, Including Two Plasmids.

Escherichia coliO157:H7 EDL933, isolated in 1982 in the United States, was the first enterohemorrhagicE. coli(EHEC) strain sequenced. Unfortunately, European labs can no longer receive the original strain. We checked three European EDL933 derivatives and found major genetic deviations (deletions, inversions) in two strains. All EDL933 strains contain the cryptic EHEC-plasmid, not reported before.

Massive horizontal gene transfer, strictly vertical inheritance and ancient duplications differentially shape the evolution of Bacillus cereus enterotoxin operons hbl, cytK and nhe.

BACKGROUND: Bacillus cereus sensu lato comprises eight closely related species including the human pathogens Bacillus anthracis and Bacillus cereus. Within B. cereus sensu lato, chromosomally and plasmid-encoded toxins exist. While plasmid-mediated horizontal gene transfer of the emetic toxin, anthrax and insecticidal toxins is known, evolution of enterotoxin genes within the group has not been studied. RESULTS: We report draft genome assemblies of 25 strains, a phylogenetic network of 142 strains based on ANI derived from genome sequences and a phylogeny based on whole-genome SNP analysis. The data clearly support subdivision of B. cereus sensu lato into seven phylogenetic groups. While group I, V and VII represent B. pseudomycoides, B. toyonensis and B. cytotoxicus, which are distinguishable at species level (ANI border ≥ 96 %), strains ascribed to the other five species do not match phylogenic groups. The chromosomal enterotoxin operons nheABC and hblCDAB are abundant within B. cereus both isolated from infections and from the environment. While the duplicated hbl variant hbl a is present in 22 % of all strains investigated, duplication of nheABC is extremely rare (0.02 %) and appears to be phylogenetically unstable. Distribution of toxin genes was matched to a master tree based on seven concatenated housekeeping genes, which depicts species relationships in B. cereus sensu lato as accurately as whole-genome comparisons. Comparison to the phylogeny of enterotoxin genes uncovered ample evidence for horizontal transfer of hbl, cytK and plcR, as well as frequent deletion of both toxins and duplication of hbl. No evidence for nhe deletion was found and stable horizontal transfer of nhe is rare. Therefore, evolution of B. cereus enterotoxin operons is shaped unexpectedly different for yet unknown reasons. CONCLUSIONS: Frequent exchange of the pathogenicity factors hbl, cytK and plcR in B. cereus sensu lato appears to be an important mechanism of B. cereus virulence evolution, including so-called probiotic or non-pathogenic species, which might have consequences for risk assessment procedures. In contrast, exclusively vertical inheritance of nhe was observed, and since nhe-negative strains appear to be extremely rare, we suggest that fitness loss may be associated with deletion or horizontal transfer of the nhe operon.

Draft Genome Sequence of Bacillus cytotoxicus CVUAS 2833, a Very Close Relative to Type Strain NVH 391-98 Isolated from a Different Location.

We report the draft genome sequence of Bacillus cytotoxicus CVUAS 2833, isolated from potato puree in Germany (2007), which is-despite its clearly different source-very similar to the type strain B. cytotoxicus NVH 391-98 isolated in France (average nucleotide identity, 99.5%).